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ABSTRACT 

Oscillations of the Sun and solar-like stars are believed to be excited stochastically by convection near 
the stellar surface. Theoretical modeling predicts that the resulting amplitude increases rapidly with 
the luminosity of the star. Thus one might expect oscillations of substantial amplitudes in red giants 
with high luminosities and vigorous convection. Here we present evidence that such oscillations may 
in fact have been detected in the so-called semiregular variables, extensive observations of which have 
been made by amateur astronomers in the American Association for Variable Star Observers (AAVSO). 
This may offer a new opportunity for studying the physical processes that give rise to the oscillations, 
possibly leading to further information about the properties of convection in these stars. 

Subject headings: convection - supergiants - stars: variables: other - stars: oscillations 



1. INTRODUCTION 

There is strong evidence that the observed solar oscilla- 
tions arc intrinsically stable but excited stochastically by 
the turbulent convection in the near-surface layers. Stabil- 
ity calculations that take the interaction with convection 
into account generally find linear stability (e.g. Balmforth 
1992a); also, it would be difficult to account for the ob- 
served very low amplitudes if the modes had been unstable 
(Kumar & Goldreich 1989). On the other hand, calcu- 
lations based on a mixing-length formulation (Balmforth 
1992b) or detailed hydrodynamical simulations of convec- 
tion (e.g. Stein & Nordlund 1998, 2001) have obtained 
energy input from stochastic convective excitation that ap- 
pears consistent with the observed oscillation amplitudes 
and damping rates. Also, the observed distribution of 
mode amplitudes (e.g. Chaplin et al. 1997) is essentially 
consistent with theoretical expectations assuming this ex- 
citation mechanism (Chang & Gough 1998). 

On this basis one may try to predict amplitudes of sim- 
ilar oscillations in other stars. Early simple estimates by 
Christensen-Dalsgaard & Frandsen (1983) indicated that 
the amplitudes increase with effective temperature along 
the main sequence, and also increase rapidly with increas- 
ing luminosity. These estimates have recently been con- 
firmed by more detailed calculations, although still based 
on a mixing- length formalism, by Houdek et al. (1999). 
In the present context it is particularly noteworthy that 
Christensen-Dalsgaard & Frandsen predicted quite sub- 
stantial amplitudes for stars on the red-giant branch. 

The predicted amplitudes for main-sequence stars are 
very low; a few parts per million in intensity or of or- 
der lms" 1 in radial velocity. Thus it is hardly surprising 
that little evidence is available for solar-like oscillations in 
stars near the main sequence. Brown et al. (1991) found 
indications of solar-like power in Procyon. This was re- 
cently confirmed by Martic et al. (1999) and Barban et al. 
(1999) who, however, noted that the inferred amplitude 
was smaller by about a factor of 3 than the predictions 



of Houdek et al. (1999). Kjeldsen et al. (1995) detected 
oscillations in the equivalent widths of Balmer lines in the 
subgiant 77 Boo, at amplitudes which were not inconsis- 
tent with the theoretical expectations; however, Brown 
et al. (1997) were unable to confirm these oscillations 
in Doppler-velocity observations. Gilliland et al. (1993) 
made an extensive search for solar-like oscillations in stars 
in the open cluster M67, but with no positive detections; in 
several cases they obtained upper limits to the amplitudes 
substantially below the theoretical predictions. For the so- 
lar near- twin a Cen A Kjeldsen et al. (1999) failed to find 
definite oscillations, although the upper limit was again 
consistent with expectations. However, Schou & Buzasi 
(2001) found clear evidence for oscillations in photomet- 
ric observations of a Cen A from the WIRE satellite, with 
amplitudes slightly higher than solar, as expected. Finally, 
Bedding et al. (2001) very recently made a definite detec- 
tion of solar-like oscillations in the sub-giant f3 Hyi, while 
Bouchy & Carrier (2001) obtained a strikingly clear signal 
of oscillations in a Cen A, in both cases with amplitudes 
close to the expected. Evidently, the observational situa- 
tion for stars near the main sequence is improving, while 
the theoretical understanding is still rather uncertain. 

For giant stars, larger amplitudes are expected. Indeed, 
Arcturus (of spectral type K1.5III) has been found to show 
oscillations that are plausibly of solar-like nature (Smith, 
McMillan & Merline 1987; Innis et al. 1988; Merline 1998). 
Similar oscillations have also been found in luminosity ob- 
servations of K giants in 47 Tuc (Edmonds & Gilliland 
1996), and in radial- velocity observations of a sample of 
K and M giants and supergiants (Cummings et al. 1998). 
Finally, Buzasi et al. (2000) found possible evidence for 
solar-like oscillations in the moderate red giant a UMa. 

For even more luminous red giants quite substantial 
amplitudes, possibly reaching several tenths of a magni- 
tude, are expected. Studies of these stars are complicated, 
however, by the fact that periods of the order of sev- 
eral months or even years are expected, requiring very ex- 
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tended observations which are generally not available from 
professional investigations. It is therefore extremely fortu- 
nate that the American Association of Variable Star Ob- 
servers (AAVSO) has a very extensive database on large- 
amplitude, long-period pulsators. In this letter, we analyze 
a sample of such stars to argue that they do indeed show 
evidence for stochastically excited, solar-like oscillations. 

2. OBSERVATIONS OF SEMIREGULAR VARIABLES 

The observations were obtained from the AAVSO 
archives, in most cases extending over 30 - 40 years; al- 
though more extensive data are available, these are so far 
the longest sequences that have been validated and an- 
alyzed in detail. The CLEANest algorithm (Foster 1995) 
was used, involving least-squares fitting of the data to sinu- 
soids taking into account the possible presence of multiple 
periods. The data were analyzed over windows of typical 
lengths of a few thousand days, hence providing measures 
of the periods P(t) and amplitudes A(t) as functions of 
time t. From this we can in particular determine the mean 
(A) and the root-mean-square standard deviation a (A) of 
the amplitudes. 

Some results are illustrated in Fig. 1. As noted already 
by Mattei et al. (1997) there is a clear division in this dia- 
gram, between the Mira variables with large amplitudes 
and small amplitude scatter and the semiregular (SR) 
variables with relatively low amplitudes and large scatter. 
Also, it is striking that there appears to be a strong corre- 
lation, in the latter group, between A and <j(A), whereas 
no such correlation is found for the Miras. This correla- 
tion is the main piece of evidence which we consider in the 
following. 
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damping time, the energy distribution is exponential, lead- 
ing to the following distribution for the amplitude A: 



(1) 



where the average pulsation energy is proportional to 
A 2 . The observed amplitudes of solar oscillation sat- 
isfy this distribution quite accurately (e.g. Chaplin et al. 
1997). According to Eq. (1) the average (A) and stan- 
dard deviation <r(A) of the amplitude are determined by 
(A) = {V^/2)A, a 2 (A) = (A 2 ) - (A) 2 = (1 - ir/4)A 2 , and 
hence 



a(A) 



1/2 

1 I (A) ~ 0.52(A) . (2) 



To illustrate this behavior, we have computed a set of 
artificial data simulating a stochastically excited oscilla- 
tor with a period of 82 days. The simulation is in addition 
characterized by the sampling-time interval At, the damp- 
ing time Td and the average excitation 5 A at each sampling 
time. The signal is written as A\{t) sm{wt)+A2{t) cos(wi), 
where u) is the angular frequency corresponding to the pe- 
riod. At each sampling time the amplitudes A\ and A2 
are reset, according to 



Ai(t + At) = exp(-At/T d )Ai(t) + 5 A Mi 



(3) 



where M% is a normally distributed random variable with 
zero mean and unit variance. The 'data' cover a total pe- 
riod of around 1600 years. They were analyzed in segments 
of 1 year, to determine the total power V% for each segment 
which was used as an estimate of the mode energy; the en- 
ergy, normalized to its average, was then determined as 
Ei/(E) = Vi/(V), where the average is over the segments. 
The mean and standard deviation of the amplitudes, esti- 
mated by y/Vi, were similarly evaluated. 

Figure 2 shows an example of a simulated timestring, 
together with the distribution of inferred mode energies, 
binned and suitably normalized (cf. Chang & Cough 1998) 
to make it directly comparable with the exponential dis- 
tribution exp(—E/(E)) which is also shown. It is evident 
that the data do indeed approximately satisfy the expo- 
nential distribution of energy. For this sample, we find 
a(A)/(A) = 0.49, in reasonable agreement with Eq. (2). 



Fig. 1. — Properties of long-period variables, from AAVSO ob- 
servations analyzed with the CLEANest algorithm. The abscissa 
is amplitude A in visual magnitude, and the ordinate is the root- 
mcan-square scatter in the amplitude. Crosses and diamonds show 
stars that have been classified as Mira variables and semiregular 
variables, respectively. Adapted from Mattei et al. (1997). 

3. PROPERTIES OF STOCHASTICALLY EXCITED 
OSCILLATORS 

The statistical properties of stochastically excited oscil- 
lations were discussed, for example, by Kumar, Franklin 
& Goldreich (1988) and Chang & Gough (1998). If the 
modes are observed over a time short compared with the 



4. RESULTS AND DISCUSSION 

According to Eq. (2) there is a direct relation between 
the value and scatter of the amplitude for stochastically ex- 
cited oscillators. This is clearly reminiscent of the correla- 
tion observed in Fig. 1 amongst the SR variables. To make 
a more detailed comparison, we show results for these vari- 
ables again in Fig. 3, together with the relation (2). It is 
evident that the observed and predicted behavior are in- 
deed quite similar. Thus it is plausible that the observed 
oscillations are excited in this manner. 
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Fig. 2. — Artificial time series for an oscillation with a period of 
82 days, a damping time of 60 years and a sampling-time interval 
At of 20 days, computed according to Eq. (3). The top panel shows 
the computed time series which, as indicated, covers about 1600 
years. In the bottom panel the points show the binned, normalized 
distribution of mode energy, in units of the mean energy; the line 
corresponds to the expected exponential distribution in Eq. (1) (for 
further details, see Chang & Gough 1998). 

In the case of the Sun, and the carefully observed K 
giant Arcturus, a substantial number of modes are ob- 
served to be excited, as would indeed be expected given the 
broad-band nature of the stochastic forcing. Christensen- 
Dalsgaard & Frandsen (1983) noted that the number of 
available modes decreases in very luminous giant stars 
where the acoustical cut-off frequency, which likely pro- 
vides an upper limit to the observable spectrum, ap- 
proaches the frequency of the fundamental radial mode. 
Even so, it is encouraging that most of the SR variables 
appear to be at least double-periodic (Mattei et al. 1997). 

Further analyses are required to reach a definite conclu- 
sion on the nature of these stars. Additional scatter in 
the inferred amplitudes might arise from beating between 
unresolved frequencies in the spectra of the stars. On the 
other hand, it was noted by Kumar et al. (1988) that the 
amplitude distribution becomes more peaked, reducing the 
scatter, if the observing time is not short compared with 
the (in this case entirely unknown) lifetime of the modes. 
An additional effect may arise because of the large ampli- 
tudes of the oscillations. Indeed, it is likely that for the SR 
variables, as for the Mira variables, the amplitude in the 
visible region of the spectrum is strongly enhanced by at- 



mospheric effects. By comparing the AAVSO amplitudes 
with the bolometric amplitudes obtained from Whitelock, 
Marang & Feast (2000), for the stars common between 
the two samples, we find that the visual amplitudes are 
on average larger by a factor 5-6 than the bolometric 
amplitudes. This may cause a nonlinear relation between 
the underlying pulsations and the observed response, pos- 
sibly involving saturation at the highest amplitudes. To 
model this we have carried out additional analyses of the 
simulated timestring shown in Fig. 2a, after filtering with 
a suitably scaled sinh -1 response. It was found that such 
saturation reduced the scatter quite substantially: even if 
the maximum amplitude was reduced by only 17 per cent, 
relative to the linear response, a (A) /(A) was reduced from 
the value of 0.49 found in the linear case to 0.47. There is 
indeed some hint that the observed values are below the 
linear trend at higher amplitudes, in accordance with this 
mechanism. 
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Fig. 3. — Root-mean-square scatter in the amplitude against 
amplitude, for the low-amplitude sample of AAVSO stars (see also 
Fig. 1). Crosses and diamonds show stars that have been classified 
as Mira variables and semiregular variables, respectively. The solid 
line shows the theoretical relation for stochastically excited oscilla- 
tors, from Eq. (2). 

The amplitudes found for the semiregular variables 
range as high as around 1.5 magnitude in the visual (cf. 
Fig. 3). Assuming the scaling between visual and bolo- 
metric amplitudes mentioned above, this corresponds to 
an amplitude in relative luminosity variations 8L/L of 
around 0.3, roughly five orders of magnitude higher than 
observed in the Sun. No estimates have so far been made 
of amplitudes of specific stars in this group. Based on the 
calculations by Christensen-Dalsgaard & Frandsen (1983), 
Kjcldsen & Bedding (1995) estimated that SL/L scales 
as LM _1 T~^ , where L, M and T c g are luminosity, mass 
and effective temperature of the star. The validity of ex- 
trapolating this relation to the parameter range of the 
semiregular variables is perhaps somewhat doubtful; also, 
the parameters of these stars, particularly their mass, are 
uncertain. However, the amplitudes observed are at least 
not obviously inconsistent with the estimates based on the 
scaling relation. 
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5. CONCLUSIONS 

At present we have secure observational information 
about very few stochastically excited pulsating stars, and 
hence additional information about this kind of pulsation, 
for a broad range of stars, would be extremely valuable, 
as a test for our models of the excitation. With a better 
physical understanding of the excitation mechanism the 
properties of the oscillations may also provide interesting 
diagnostics of the properties of convection, as a supple- 
ment to direct observations of velocity and intensity effects 
of surface convection. Secure identification of stochasti- 
cally excited pulsations in the scmircgular variables would 
therefore evidently be of substantial interest. 

The present study is of a statistical nature. To confirm 
our suggestion a more detailed analysis of individual stars, 
including careful modeling of the stellar structure and the 
properties of the oscillations, will be required, furthermore 
using the frequencies to constrain the overall properties 
of the stars. Also, further observational data should be 
sought, making full use of the AAVSO data archive, and 
of other similar archives of long-term observations. In- 
deed, the century-long data series available for a number 
of stars in the AAVSO archive are now being validated and 
quality controlled; these full datasets will then be submit- 
ted to the same type of analysis as reported here. To 



obtain additional observations of these long-period vari- 
ables will clearly be a long-term project, although one 
definitely worth undertaking. However, it is also of evi- 
dent interest to study similar variations in stars further 
down the red-giant branch, such as already observed in 
K giants, of shorter periods. Here the predicted ampli- 
tudes are also much smaller, precluding visual observation 
of the variations. However, with the present trend towards 
increasing sophistication of amateur equipment, including 
the extensive use of CCD detectors, there is little doubt 
that amateur astronomers will be able to continue and 
extend their very valuable efforts towards providing long- 
term observations of a nature that are difficult to obtain 
in a professional setting. 
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